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WAVE EQUATION ON ONE-DIMENSIONAL FRACTALS WITH SPECTRAL 
DECIMATION AND THE COMPLEX DYNAMICS OF POLYNOMIALS 

ULYSSES ANDREWS, GRIGORY BONIK, JOE R CHEN, RICHARD W. MARTIN, 

AND ALEXANDER TEPLYAEV 


Abstract. We study the wave equation on one-dimensional self-similar fractal structures that can 
be analyzed by the spectral decimation method. We develop efficient numerical approximation tech¬ 
niques and also provide uniform estimates obtained by analytical methods. 
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1. Introduction 

The purpose of this paper is to study, both analytically and numerically, the wave equation on 
the unit interval endowed with a self-similar fractal measure. Previous studies of wave equation on 
fractals, including numerical approximations, were published in |l20l|24l|25l|2S|33l|59]|. All these 
works have some, although not direct, relation to the classical paper Il55]l . but are more directly 
related to the fractal Fourier analysis, see Il57ll58l . Our computational methods mostly come from 
the theoretical papers |l9l[l0llSll62]] that develop so-called spectral decimation method in the form 
applicable for to numerical analysis. 

In general, there is a large literature dealing with analysis and probability on fractals in math¬ 
ematical terms, such as |1 [HI [12 [13 [II [M EHl [32l [M |35l E Ell 0111^ llQl |45l |5l 123 |60l 
and references therein], and also extensive mathematical physics literature, including [[TIEIISIEIEII 
123 ESI [23 EH 113 mi- Of particular interest are the works studying the appearance of fractals in 
quantum gravity, including 01 HI El [13 [El [13 [23 E31131111 [121 [50[j . 
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We consider a situation in which a good enough (fractal) Laplacian A is defined on L‘^{K, ji), 
where a compact set K (the unit interval in our case) equipped with a (fractal) Borel measure /i. 
This Laplacian A is a point-wise limit or as the generator of a Kigami’s resistance form (see Propo¬ 
sition |2.3[ ), and one can extend some of the classical numerical techniques to approximate some 
(intrinsically smooth) solutions of the wave equation initial value problem 

( dtfU = —Am on K X [0,T], 

(1.1) < m(-, 0) = (p on K, 

y dtu(-, 0) = Ip on K. 

As is well known, if the spectrum of the Laplacian is discrete, then the solution of the wave equation 
can be represented in terms of L^(;u)-eigensolutions {A^, fk}^=o of the Laplacian A, with Aq < 
Ai < A 2 < • • • t oo and Afk = Xkfk- Writing 

00 00 

(1-2) ^ = XI ^ = X 

k=0 k=0 

where = (0, fk)L^ and 13k = {ip, fk) l^, one finds that u admits the series representation 


(1.3) 


u{x,t) = y^^akfk{x)cos 

k=0 


00 

k—k ■ 

IX —O/jjiin 


^^/fc(a:)sin 
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where /Cmin := min{j G N U {0} : Xj > 0}. It is known that the series point-wise converges 
poorly and the numerical approximations are very unstable unless the smoothness of solutions can 
be controlled. 

In our setup, K = I and /r is the fractal measure defined in Section For simplicity we assume 
that the initial velocity xp = Q, and so the solution to (|L1|) is 


(1.4) 


u{x, f) = X ^kfk{x) COS U\A 


k=0 


If we theoretically assume that cp is given by a (^-impulse at point 0, 0 = then we have that 
:= fk{x)So{x) ii{dx) = /fc(0). Note that (5o(a:) is not a function by the unit atomic measure 
at zero, and so the integral in this definition of is to be understood as a formal expression, as 
in the theory of distributions (for the classical version, see ll56ll . and for the fractal version, see 
II 52 II L This approach on a fractal space does not allow an accurate numerical approximation of the 
solutions. 

Therefore we concentrate on a situation where the initial condition is highly localize function, but 
is smooth in intrinsic sense, and we can show that the approximating series converges uniformly. 
This is an illustration of the general principle of Stricharz ||58]|: Laplacians on fractals with spectral 
gaps have nicer Fourier series. However, the abstract result [|5^ does not include the estimate of 
the remainder which we obtain in our work. 

Numerically, we can only compute the eigensolutions of the fractal Laplacian up to a finite level, 
so in practice we solve the “approximate” wave equation 

^tt^n AjiXiji on Vji X [ 0 , T], 

= onK, 
fi) fi 


(1.5) 
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where 0'kfn,k is the approximate (^-function built up from the first |14(,| = (3”° + 1) 

k=0 

eigenfunctions of A„ (with Anfn,k = ^n,kfn,k), and ak := «no,fc > 0 are the coefficients found in 
Section |33| Throughout the section no will be fixed, and we will not mention uq explicitly unless 
the context demands it. The solution to ( |1.5| ) has the series representation 

(1.6) Un{x,t) = y^^akfn,k{x)cos forallx e K andf e [0,T]. 

k=0 

For each t, we harmonically extend the function x Un{x, t) from 14 to /. This procedure allows 
us to compare Un{x, t) with 

3’"o 

(1.7) u{x, t) = oikfk{x) cos for all a: G / and t G [0, T], 

k=0 

the solution of the wave equation on (/, n) whose initial condition is the truncated series represen¬ 
tation of the (5-impulse. We note that u is differentiable in t and continuous in x. However it is 
highly localized function at f = 0, and therefore it mimics wave propagation from a delta function 
initial values. 

Our paper is organized as follows. Section contains the construction of the unit interval as a 
p.c.f. fractal, definition of the Dirichlet energy form, the definition of the corresponding Laplacian 
and its associated eigenvalues. In Section we use spectral decimation to construct the eigen¬ 
functions of the discrete Laplacian and prove that their limit is continuous. The section concludes 
with the spectral decomposition of the delta function. Section contains various technical es¬ 
timates needed to show the convergence of solutions of the wave equation. In Section we give 
theoretical bounds on the approximations to the wave equations solutions and convergence informa¬ 
tion. Section [^contains the numerical computation of the wave equation solutions, their associated 
eigenfunctions, and the Fourier approximations for the delta function. 


Remark 1.1. Theoretically, the infinite propagation speed for wave equation solutions was estab¬ 
lished in on some p.c.f. fractals with heat kernel estimates 


( 1 . 8 ) 


Cl 


V(x, 


exp -C2 


d(x,yy \ 

t ) 




< p{t,x,y) 


< 


C3 




exp — C4 


dix,y)>^ \ 




for positive constants Cl, C2, C3, C4, a:, 1 / E I,t E (0,1], where /5 = 2/ds and V{x,r) = ii{Br{x)). 
Kigami in ll^ obtained such estimates in a situation which resembles, but is technically different, 
from ours. We conjecture that an analogue ( |1.8| ) holds in our situation, but proving this would lie 
outside of the scope of our paper. 


Acknowledgement. The authors are very grateful to Daniel Kelleher, Hugo Panzo and Antoni Br- 
zoska for many helpful discussions, and to Luke Rogers for explaining the eigenfunction estimates 
based on his paper llAfll . A.T. also thanks Sze-Man Ngai and Alexander Grigor’yan for very valuable 
advice. The authors thank anonymous referees for corrections and a substantial list of constructive 
suggestions leading to improvements in the first version of our paper, and for the suggestion to 
include the infinite wave propagation speed Remark [Llj 
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2. Eigenvalues of the fractal Laplacian on an interval 


In this section we define a particular self-similar structure on the unit interval. In this way, it 
can be seen as a p.c.f fractal (see IfTTl [T 6 l |33 [37l |Sl |63l)- In these papers the reader can find 
these definitions and an exposition of the general theory of Dirichlet forms on fractals, as well as 
further references on the subject. Herein we will use three contractions for simplicity. However, 
one could perform the same construction using any number of contractions in order to obtain a 
fractal Laplacian on the unit interval. 

To define the standard Laplacian, we can use three contractions Fi, F 2 , F 3 : M ^ M Fj{x) = 
\x -\- ^pj with respective fixed points pi = 0, p 2 = Pa = 1- Then the interval /=[0,1] is a unique 
compact set such that I = IJ Fj{I). The discrete approximations to I are defined inductively by 

i=1.2,3 

K = U FjiVn-i) = where Vq = dl = {0,1} is the boundary of /. For x,y eVn we 

J ~ 1 5 2 , 3 

write y x if \x — y\ = S”**. Then the standard discrete Dirichlet (energy) form on I 4 is 

£n(/,/) = 3” {f{y)-f{x)f, 

x,yeVn 

yr^x 

and the standard Dirichlet (energy) form on / is £(/, /) = lim £n(/, /) if this limit exists. We 

n—^-OO 

call a function h harmonic if it minimi z es the energy subject to the constraint of the given boundary 
values. Then we have that £„+i(/, /) ^ £«(/, /) for any function /, and £„+i(/j, h) = 8,n{h, h) = 
E{h, h) for a harmonic h. A function h is harmonic if and only if it is linear. If / is continuously 
differentiable then 

£(/,/)= f\f{x)\^dx. 

Jo 

The domain T of this standard Dirichlet (energy) form £ on / coincides with the usual Sobolev 
space H^[0, 1]. Moreover £ on / is self-similar in the sense that 


£(/,/) = 3 

i=l,2,3 

The corresponding standard discrete Laplacians on I 4 

^nf{x) = ^ if{x) - f{y)), X e Vn \ Vo , 

y^y-n 

yr^x 


and the (renormalized) Laplacian on / is 

Mix) = lim TMfix) = -I fix) 

n^oo Z 


for any twice differentiable function. In our convention the Laplacian is a nonnegative operator. 
For any twice differentiable function /, the Gauss-Green (integration by parts) formula applies 

8,{fJ) = ^ t f^fdx + ff\ 

Jo 0 

We can modify the above construction with the introduction of the parameter p, where 0 < p < 1, 
and write g = 1 — p. Later we will show that these parameters give the transition probabilities of a 
random walk on the unit interval. Now, we define contraction factors (or resistance weights) 


P 

1 +p 


9 

1 +p’ 


( 2 . 1 ) 


ri = rg = 


and r 2 = 





and measure weights 
( 2 . 2 ) 
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^ A P 

nil = m 3 =- and m 2 = -. 

1 + q l + q 

Note that in general the choices of resistance and measure weights are essentially free, up to constant 
multiples, according to Kigami’s theory of Harmonic calculus on p.c.f. self-similar sets (361 [37]], 
but we make a unique choice that leads to a manageable spectral analysis, as explained in (9115311601 
I6lll62ll6l. We do not give a complete explanation here because it would require too much space. 
In short, the spectral decimation requires a symmetry mi = m 3 . Moreover, the spectral decimation 
also requires that the resistance weights are, up to a constant, reciprocals of the measure weights, 
and 

(2.3) nil + m 2 + m 3 = ri + r 2 + n = 1. 

Thus, our system essentially has one independent parameter, which we denote p and express every¬ 
thing else in terms of this parameter. 

We may now define the three contractions: Fi, F 2 , F 3 : M —> M with respective fixed points 
Pi = 0 ,P 2 = = 1 in terms of resistances which depend on our parameter p 

(2.4) Fj{x) = TjX + (1 - rj)pj. 

Then the interval /=[0,1] is the unique compact set such that 

(2.5) 1 = 

i=L2,3 

The discrete approximations to / are defined inductively by 

( 2 . 6 ) K = U i"i(K-i), 

i=l,2.3 

where Vo = 5/ = { 0 ,l}is the boundary of /. 

The following definitions and results come directly from the more general theory in (miT6ll36l 
[ 37 I . so we omit the proofs. 

Definition 2.1. The discrete Dirichlet (energy) form on I 4 is defined inductively 

(2.7) £n(/,/)= 5] ^En-l{foFj). 

i=1.2,3 

with £o(/) /) = (/(I) “ /(0))^> nnd the Dirichlet (energy) form on I is 

(2.8) £(/,/)= lim £„(/,/)= /" \f{x)\^dx 

n^oo Jq 

The domain 5" of £ consists of continuous functions for which the limit is finite, and coincides with 
the usual Sobolev space H^[0, 1]. 

The existence of this limit is justified by the next proposition. 

Proposition 2.2. We have that £^+i(/, /) ^ £n(/) f) for any function f, and 

(2.9) E^+i{h,h) = Enih,h) = E{h,h) 
for a harmonic function h. 
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Proposition 2.3. The Dirichlet (energy) form E on I is local and regular, and is self-similar in the 
sense that 


(2.10) £(/,/)= 5] ^EifoFjJoFj). 

i=l, 2,3 


The domain of E, see Definition \2.1\ is dense in the space of continuous functions on I. 

The p-Laplacian A^, satisfying the following Gauss-Green (integration by parts) formula 

(2.11) UfJ) = C f fAJdn + ff'll, 

Jo 

where p, is a unique probability self-similar measure with weights mi, m 2 , m 3 , that is 


( 2 . 12 ) 

can be defined by 
(2.13) 


l^= mjpoFj. 

i=l, 2,3 

A^,/(x) = lim (l+^)’'A„/(a:), 


where the discrete Laplacians 


(2.14) E^nf{Xk) 


f{xk) -pf{xk-i) - qf{xk+i) 

or 

f{xk) - qfixk-i) - pf{xk+i) 


are defined as the generators of the nearest neighbor random walks on I 4 transition proba¬ 
bilities p and q assigned according to the weights of the corresponding intervals. The domain of 
the corresponding continuous Laplacian A^, defined to be the set of all continuous function f for 
which the limit ( |2.13| ) exists and is continuous, is dense in the space of continuous functions on I. 


Note that by definition p = — and q = —. The transition probabilities p and q can be 

^ mi+m2 ^ mi+m2 ^ jt 'i 

assigned inductively as shown on Figure |2Tj 


Proposition 2.4 (Self-similarity of the Laplacian). 


(2.15) 


A^(moF^) 



-kl 

(A^u) O Fyj. 


•- 






-• 

1 

mi 


m2 


m3 

1 

1 


■g p 


p'g 


1 


1 q p p q q p q P P q P q q p p q 1 


Figure 2 . 1 . Random walks corresponding to the discrete Laplacians A„. 

The above construction of the standard Laplacian and the associated Dirichlet form on / corre¬ 
sponds to the case p = 1. In the p 7 ^ 1 case, a change of variables can either turn the Dirichlet form 
into the standard one, or turn the p-measure into Lebesgue measure, but not both at the same time. 
For this reason, different values of p give different p-Laplacians even up to a change of variable. 
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We can apply the classical result of Kigami and Lapidus to show that both the Dirichlet and 
the Neumann Laplacians satisfy the spectral asymptotics 

(2.16) 0 < liminf - 7 ^^ < limsup < 00 , 

A^oo \ds/^ A^oo 


where as before p{\) is the eigenvalue counting function, and the spectral dimension is 
‘ ’ ''•-log ( 1 + 1 )-'’ 


where the inequality is strict if and only if p 7 ^ g. 

In the lemma below, (t(A„) is the spectrum n of the level n Laplacian A„. 


Lemma 2.5. If z 1 ± p, then R{z) G cr(A^) if and only if z G (j(A„+i), with the same multi¬ 
plicities. Here 

z(z^ — 3z -\-2 + 
pq 

Moreover, the Neumann discrete Laplacians have simple spectrum with cr(Ao) = {0, 2} and 

n—1 

(2.19) a{An) = { 0 , 2 } (J R-^{1 ± q} 

m=0 


(2.18) 


R(z) 



for all n > 0. In particular, for alln > 0 we have 0,1 ± g, 2 G cr( A„). Also, for alln > 0 we have 
1 ± p G cr(A„) if and only ifp = g. 


Proof In this case, according to IIMl Lemma 3.4], Il45l (3.2)], we have that R{z) 
po and Pi solve the matrix equation 

(2.20) S-zIo- X(Q - zIi)-^X = po(z)Ho - pi{z)h. 

with S = Iq = Ii — 12x2, X — —g/ 2 x 2 i X = —12x2, 

1 -p 


( 2 . 21 ) 

and 


Q = 


-p 1 


( 2 . 22 ) 





, where 


Solving this we obtain 

(2.23) 
and 

(2.24) 


^o(z) 


Pi(^) = 


^ PQ 

— 2z + 1 — p2 

z(z^ — 3z + 2 + pq) 

— 2z + 1 — p2 


Then we use the abstract spectral self-similarity results (see ll^ |45]|) to find that (t(A„_|_i) = 
i?“^{(T(A„)}. Note that 0 and 2 are fixed points of R{z). The preimages of 0 are 0, 1 -I- p and 
1 + g. The preimages of 2 are 2,1 — p and 1 — g. If p 7 ^ g then 1 ± p are not eigenvalues because 
they are poles of po{z) (see ll6ll|45]|). □ 
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Figure 2.2. Sketch of the cubic polynomial R{z) associated with the fractal Lapla- 
cians on the interval. 




Figure 2.3. Eigenvalues of the first three fractal levels arranged in a rooted tree and 
numbered in increasing order. 


Remark 2.6. In Figure we give a sketch that describes the complex dynamics of the family of 
cubic polynomials associated with the fractal Laplacians on the interval (see IMIIIS]]). The curved 
dotted line corresponds to the case when p = \ and the Julia set is the interval [0,2]. For any other 
value of p (0 < p < 1, p 7 ^ |), the graph of the polynomial R{z) behaves like the shown solid 
curved line. It is easy to see that then the Julia set of R{z) is a Cantor set of Lebesgue measure 
zero. Note that the transformation p i-A 1 — p does not change the polynomial R{z), although the 
Laplacians are different. 


3. Spectral Decimation and Eigenfunction Approximations in the limit 

Thus far we have described the spectral decimation which allows us to characterize the eigenval¬ 
ues of the fractal Laplacian (Lemma [A5] ). We now turn to the eigenfunctions. 

3.1. Eigenfunction extension. In this subsection we demonstrate how to extend an eigenfunction 
to an eigenfunction fn+i,* using spectral decimation. 

To fix notation, let xq < yo < yi < xi be four consecutive vertices in I 4+1 with xq, G I 4 and 
yo,yi e I4+I \ Kx. Given an eigenfunction of A„ with eigenvalue A^,*, we define its extension 
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0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 


0 0.2 0.4 0.6 0.8 1 



Figure 3 . 4 . Fourier approximations for delta function, p = |. Left to right: no = 2,3,4. 


fn+i,* to Vn +1 according to the formulas 


(3.1) 


(3.2) 


4n+l,*(?/o) 


/n+l,*(?/l) 


9(1 - z)fn,4xo)+pqfn,*{Xi) 
{1 - p - z){l + p - z) 

9(1 - z)fn,*{Xi)+pqfn,*{Xo) 


{1 - p — z){l + p - z) 

Here we assume z 7 ^ 1 ± p. The claim is that fn+i,* is an eigenfunction of A^+i with eigenvalue 
z = R~^{Xn ,*), where R is the cubic polynomial which appeared in Lemma 2.5 As explained in 


the proof of 2.5, the preimage i?“^([0, 2 ]) has three branches, so each eigenvalue on level n 
generates three new eigenvalues 4 , on level (n + 1). The only exceptions are the eigenvalues 0 
and 2, each of which generates two new eigenvalues because 1 ± p are forbidden (see Figure [23] ). 
This means that each eigenfunction extends to either two or three eigenfunctions at the next level. 

Theorem 3.1 (Eigenfunction extension). Suppose : I 4 ^ E A an eigenfunction of An with 
eigenvalue Xn,*. Let fn+i,* ■ 14+i M.be an extension of fn,* to Vn+i defined via (3.1) and (3.2), 
with z 1 A p. If 

z^ — 3z^ + (2 + pq)z 


(3.3) 


Xn,* ■^(^) 


P9 


then fn+i,* is an eigenfunction of An+i with eigenvalue z = R ^(A„^*). 

Proof We break the proof into two parts. Given fn,*, we first show that the following are equivalent 
for an extension fn+i,* of fn,*'- 

(1) fn+i,* is defined via the extension formulas ([34]) and (3.2). 

(2) fn+i,* satisfies the eigenvalue equation An+ifn+i,* = zfn+i,* on I 4+1 \ K- 

After establishing this equivalence, we proceed to show that fn+i,* is an eigenfunction of on 
all of Vn+i, provided that ( |3.3| ) holds. 
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Xn 


y'o 


X 


y[ 


x'l 


P Q 


b a 


q p 


Figure 3.5. A diagram of two adjacent {n + l)-cells used in the proof of Theorem [J!Tj 


First we show the equivalence of (1) and (2). Assuming (2), we apply the eigenvalue equation 
^n+i/n+i,* = zfn+i,* at the points yo, yi G Vn+i \ K to obtain, by using both formulae in (2.14) 
depending on the point x (in fact, to cover these cases as well as the case x &Vq, i.e. x is a boundary 
vertex, below we use the parameters a G {p, q, 0,1}, and b = 1 — a, instead of p and q = 1 — p). 


(3.4) 

(3.5) 


(1 - z)fn+l,*iyo) = Pfn+lAyi) + Qfn+l,*{^o), 
{1- Z)fn+l,*{yi) = qfn+lA^l) +Pfn+lAyo)- 


This is a linear system of 2 equations with 2 unknowns (/n+i,*(To) and fn+i,*Ai) are known, 
fn+i,*{yo) and fn+i,>f{yi) are unknown), whieh has a unique solution. After some elementary eal- 
eulation, and using the faet that /n+i,=^|v, = fn,*, h is easy to verify that fn+iAvo) and fn+iAvi) 
are uniquely expressed in terms of fn,*{xo) and fn,*{xi) according to the extension formulas ( |3.1[ ) 
and p.2| ), which shows (1). The reverse implication (1) (2) is straightforward. 

At this point we have proved that the eigenvalue equation An+ifn+i,*iA = ^fn+iAA holds for 
X e 14+1 \ 14 . However, we have neither used the property of the eigenfunction fn,*, nor related 
to the eigenvalue A„ *. To do so we must eheck the A„+i-eigenvalue equation on I 4 . 

We introduee some additional notation. Fix an x G I 4 . Let x'q,x'-^ G I 4 be adjaeent to x on 
level n, and yQ,y'i G I 4+1 \ I 4 be adjaeent to x on level (n + 1), as shown in Figure [T5j (If 
X G 14 = (0,1}, then there is only one adjacent vertex on level n. This will be taken eare of in the 
next argument.) We also label the transition probabilities aeeording to ( |2.14| ); see also Figure [2Tj 
The parameter a can be one of {p, g, 0,1} depending on x. In particular, to take into account that 
X G 14 has only 1 adjacent vertex, we set o = 0 if x = 1 and a = 1 if x — 0. The parameter b is set 
to equal 1 — a. 


Now we show that if z lip and (3.3) holds, then A„+i/„+i^*(x) = zfn+i,*ix) for x G I 4 . 
By ( |2l4l ), 

(3.6) An-\-lfn+l,*A^ fn+l,*{x^ ttfn-\-l,*{yi) (?/o) > 

(3.7) Anfn,*ix') fn,*{A) 

Using the extension formulas ( |3.1[ ), (3.2), /„+i *(x) = fnAA^ (3.6), we find 

(^n+l A 

4(1 - AfnAA +PQfnAAy 


= {l-z)fnAx)-a 


- b 


(1 - p - z){l +p- z) 
q{l - z)fnAA+PQfnAA) 


{1 - p - z){l + p - z) 
(1 - z){l -p- z){l -\-p - z) - q{l - z) 


(1 -p - z){l +p- z) 

PQ AfnAA) + bfnAA)) 

{1 — p - z){l + p — z) 


fnAA 


(3.8) 
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Using p.7[ ) we can write 

(3.9) + bfn,*{^o) = -(A„ - l)fn,4x) = -{Xn,* “ 

This allows us to replace the second term of p.8| ), so that the entire p.8| ) equals 

{1 - z)[{l - p - z){l + p - z) - q] - pq{l - Xn,*) 


(3.10) 


(1 - p - z){l -\-p- z) 

— - Zz^ + (2 + pq)z) + pg - pg(l - A„,*) 


-fn,*{x) 




(1 - p - z){l +p - z) 

Infer that (A,i+i — z)fn+i,* = 0 on 14, and in turn on 14 + 1 ; if .s 7 ^ 1 ± p and 

z^ — + (2 + pq)z 


(3.11) 


Xn,.* 


pq 


= R{z). 


□ 


3.2. Continuity in the limit. In this subsection, we show that the eigenfunction extension al¬ 
gorithm (Theorem |3.1| ) produces a continuous eigenfunction of the fractal Laplacian in the limit 
n —7- (X), provided that one always chooses the lowest branch of the inverse map R~^ at all levels 

n > Hq. 


Lemma 3.2. Fix no, k eNU {0}. Let fno,k ■ 14o -^^be an eigenfunction 0 /A„g with eigenvalue 
Xng,k- Let {fno+i,k}i^i ^^e Sequence of An^+i-eigenfunctions extended from fno,k via successive 
applications of Theorem |3.7| where one always chooses the lowest branch of the inverse cubic 
polynomial R~^{z) [see {3.3)]. Then 


(3.12) 


limsup max \fna+i,k{.x)\ 

i—>-oo 


is bounded. 


-1 


Proof. From Lemma 2.5 we know that Xn,* = R{Xn+i,*)- Assume that the lowest branch of R 
is chosen to generate Xn+i,k = R~^{Xn,k) from Xn,k- Observe that R is concave on [0, min(p, g)]; 
therefore the graph of i? on [0, min(p, g)] lies above the secant line connecting (0,0) and (min(p, g), 2) 
(see Figure |2!2l). This implies the inequality 


(3.13) 


Xn,k R(^Xn-\-l,k) 4 


min(p, g) 


A 


nf-l,k • 


By iterating this inequality, we see that the Afold preimage A„g+j ^ — R ^{Xno,k), where the lowest 
branch of R~^ is always chosen, satisfies 


(3.14) 


Ano+ijfc A 


l' min(p, g) 

V 2 


A 


no, 


k- 


The corresponding eigenfunctions fnQ+i,k are generated via Theorem [TTj 

LetM,i_fc = maxxeVn \fn,kA)\- For each n > no and each p G I4+1 \ 14 »we use the eigenfunction 
extension algorithm p.l| ) and ( |3.2| ) to arrive at the following estimate: there exist xo,xi G I4 such 
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that 


(3.15) 

\fn-\-l,k{y)\ 

(3.16) 

< 

(3.17) 

< 


q{l - Xn+l^k)fn+lA^o) +Pqfn+lA^l) 


(1 p A72+l,fc)(l “1“ P 

9(1 - Xn-\-l^k )|/n+l,fc(a:o)| +pg|/n+i,fc(a^i)| 

(1 P '^n+l,fc)(l ~\~ P -^n+l,fc) 
-Mn^k- 


Q ^n+l,k 

In the second line we used the triangle inequality and the bound \n+i,k < which can be seen 
from p.l4| ). This then implies the estimate 

(3.18) Mn+iy < - - -M, 


n.k- 


for all n > no. Applying p.l8| ) inductively and using ( |3.14| ), we see that for all i G N, 

_ _ f min(p, q) \ ^ X^^^k 

Setting 7 j = (^ and taking the limit, we obtain 


(3.19) 


Mno+i,k < Mno,k TT - ^ - 

^ ~ ^no+i,fc 




-1 


(3.20) 


lim sup Mno+i,k < Mno,k lim sup 11 ( 1 --T.r' 


J = 1 


It remains to verify the convergence of the infinite product n^i(l “ which is equivalent 

to showing the convergence of the series log(l — 7j)“^. Observe that if we set to satisfy 
(1 — 7j)“^ = 1 + Cj, then 


(3.21) 


log(l - 7 j) ^ = log(l + ej) < 


by the inequality 1 + x < e*. Moreover, since jj = Kr^ for suitable positive constants K and 
r < |, we can always find a constant Kq such that 


(3.22) 


= 


7i 


Kr^ 


1 — 7j 1 — Kr^ 


< Kor^ 


for all sufficiently large j. Since the geometric series Kov^ converges, this implies that the series 
log(l ~ oij)~^ converges. □ 

Now we prove the continuity of the eigenfunction in the limit. 

Theorem 3.3. Let {fno+i,k}^i be the sequence of IXnQ+i-eigenfunctions extended from fno,k ns 
in Lemma 3.2 Then the limit fk ■= limj^.oo fno+i,k is uniform on I, and can be extended to a 
continuous function on I. 

Proof. The key argument is that since the eigenvalues Xn^+i^k = -R~*(^no,fc) to 0 as i —?• 00 , the 
eigenfunction extension ( |3.1| ) and ( |3.2| ) of fno+i,k to Vno+i+i can be approximated by the harmonic 
extension of fno+i,k to Vno+i+i as i 00 , uniformly on /. Since a harmonic extension on / is 
continuous in the limit, we deduce that the limit fk can also be extended to a continuous function. 
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Define, for each n > no and each k, the harmonic extension fn+i,k of fn,k to Using the 

coordinates xq, xi,yo,yi introduced before Theorem |TTj 


(3.23) 

(3.24) 


fn+l,kiyo) 

fn+l,k(jjl) 


fnA^o) + Pfn,k{^l) 
l+P 

/n,fc(^l) T Pfn,k{^0^ 
l+P 


Note that these are p.l| ) and p.2[ ) with z = 0. 

Let us now estimate \fn+i,k{yo) “ /n+i,fc(l/o)|, which equals 


(3.25) 

(3.26) 


q{l - Xn,k)fn,kM + PQfn,k{^l) _ fn,k{Xo) + pfn,k{^l) 
{1 - P - Xn,k)il + P - K,k) l+P 

+ P PQ X^j^fjikixo^ + pXji^ki^ X^k^fjikix\) 
(? '^n,/c)(l + P '^n,/c)(l 3 “ P) 


Using the trian gle ine quality and then replacing | fn,k (xo) \ and | /n,fc (xi) | by the sup Mn,k = sup^.^ I fn 
we can bound p.26| ) from above by 

-^71,11(1 + P PQ Xn^k^ ]l+n,k + pXn,k{^ Xn^k^Mn,k ^n,/c 


(3.27) 


{q - An,fc)(l +p- A„,fc) (1 + p) 


q Xjik 




Since lim^^oo Afn.fc is bounded by Lemma [T2| and A„ ^ ^ 0 as n —?• cx), the right-hand side of this 
inequality tends to 0. The same estimate holds for \fn+i,k{yi) — fn+i,k{yi)\- Since yo and yi are 
arbitrary, we conclude that \fn,k — fn,k \ converges to 0 uniformly on I. □ 


3.3. Spectral decomposition of the delta function. Let {fn,k]T=Q ^ complete set of eigen¬ 
functions of A„ with corresponding eigenvalues {Xn,k}T=o- Consider the level-n delta function 
: K ^ E defined by 


(3.28) 



if X = 0, 
if X e 14 \{ 0 }. 


By the spectral theorem, we can find a set of real numbers (or weights) {an,k}kLo such that 


gn 

(3.29) = '^an,kfn,k{^)- 

k=0 

The sequence approximates a delta function at 0 in the limit n +■ 00 . 

In order to study the wave equation in Section|^ we need estimates on the eigensolutions (/„,,*, A„^*), 
as well as information about the weights We will address the former in Section^ and the latter 
in the following proposition. 


Proposition 3.4. The weights «„+!,* can be obtained inductively from an,*- 

Proof First we fix our convention at n = 0. The two (non-£^-normalized) eigenfunctions fo,i and 
/o ,2 of Ao are 

(3.30) (/o,i(0), /o,i(l)) = (1,1) and (/o,2(0), /o, 2 (l)) = (1, -1), 

with corresponding eigenvalue 0 and 2, respectively. It is then easy to see that = lfo,i + |/o, 2 . 

i.e., ao,i = a;o ,2 = 
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For the iteration step, suppose the weights * are known at level n, and we want to determine 
the weights The idea is to write each contribution <y.n,kfn,k in terms of a linear combina¬ 

tion Y.- «n+i,fc, /n+i,fc, of the (2 or 3) eigenfunctions fn+i,kj which are extensions of fn,k given by 
Theorem 13.11 

To make this idea precise without adding too much notation, we fix n and k, and write /, A, 
and a as respective shorthands for fc, and fc. If A ^ {0,2}, then spectral decimation 

(Theorem |3.1| ) implies that / has 3 extensions /i, / 2 , and /s to Vn+i which are eigenfunctions of 
A„+i with respective eigenvalues Ai < A 2 < A 3 . We would like to find the corresponding weights 
«!, 02 , and as by imposing the following matching condition: For any four consecutive vertices 
xo<yo<yi < xi in I 4+1 with xq, xi e K and yo, yi e K+i \ K, 


(3.31) 


3 

'^Oiifiix) 

1=1 


fa/(xo), ifx = Xo, 

[a/(xi), if X G 


An explicit calculation verifies that with the weights an+i,* generated from this matching condition, 
we have 


(3.32) 


gn+l 

5 ] an+l,kfn+l,k = 4 ”^'^ 

k=0 


We now determine the weights. Observe that fi and / agree on Vn by construction. This together 
with the matching condition p.31| ) at xq (or at xi) implies that 

3 

(3.33) ^ ^ ai = a. 

i=l 


Next, using the eigenfunction extension formula p.l| ) and the matching condition at yo in ( |3.31[ ), 
we get 


(3.34) 


E 


5(1 - \) 


i=l (1 “ P “ (1 + P “ 


ai = 0 . 


Notice that there is no dependence on /. Similarly, using ( |3.2| ) and the matching condition at yi in 
( |3.31[ ), we arrive at a third relation 

3 

(3.35) y] 


pq 


rf ~ P ~ P ~ \) 


ai = 0 . 


i=l 


Equations ( |3.33[ ), ( |3.34 ), and(3.35)formalinear system of3 equations with 3 unknowns (ai, a 2 , as). 
It has the unique solution 


(3.36) 

(3.37) 

(3.38) 


(A3 — A2)(q' — Ai)(l + p — Ai) 
(A2 — Ai)(g — A3)(1+p — A3) 

(Ai — A 3 )(g — A2)(1 + p — A2) 
(A2 — Ai)(g — A3)(1 + p — A3) 

_ (1 + P ~ — A3) 

(A3 — Ai)(A3 — A2) 
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It is possible to find «3 from a and A 3 only. From (3.3) we know that i?(Aj) = A for i G {1,2,3}, 
which means that 


(3.39) (z - Ai)(z - A2)(z - A3) = pqR{z) - A. 

By diiferentiating both sides of ( |3.39| ) with respect to z, and then evaluating at z = A3, we obtain 

(3.40) (A 3 — Ai)(A 3 — A 2 ) = pqR'{Xs) = 3 A 3 — 6 A 3 + 2 + pq. 


This allows us to replace the denominator in the RHS of p.38[), which leads to 


(3.41) 


{1+p- A3)(g - A3) 

-Q/ 

3 A| — 6A3 + 2 + pq 


This proves the induction from 4 , to * in the case where the eigenvalue A ^ {0,2}. 

If A G {0, 2 }, then / has 2 eigenfunction extensions to the next level. The matching condition 
stated in p.31| ) remains the same, but degenerates to a linear system of 2 equations with 2 un k nowns. 
We omit the details. 

As a simple corollary, we now show that the weights are all nonnegative in our conve ntion . 
Recall that 00,1 = Q!o ,2 = ^- By the structure of the cubic polynomial R{z) (see Figure [2^ , 
Ai G [0, min(p, g)], A 2 G [max(p, g), min(l+p, 1 + g)], and A 3 G [max(l+p, 1 + g), 2]. So if a is 
nonnegative, it is direct to verify using ( |3.36[ ) through p.38| ) that «i, a 2 , and ^3 are all nonnegative. 
By induction we deduce that all weights are nonnegative. □ 


Usin g the aforementioned result, we now define the “approximate delta functions.” Based on 
Lemma 2.5, and the fact that the lowest branch of R~^{z) is increasing, we can deduce that the 
lowest 1 14 1 eigenvalues of A„+i are determined recursively by 

Xn+i,k = (The lowest branch of R~^){Xn,k) for 0 < < 3”^. 


(3.42) 


Given the level-n delta function Sq^\ we define its approximation at level no < n by 




(3.43) 


5i”°’”^(a:) := Y] ano,kfn,ki^) for x G I 4 . 


k=0 


In other words, we consider a truncated series of the spectral representation at level no, fixing the 
coefficients ano,k, but taking the eigenfunctions fn,k to level n. 


4. Estimates of eigenvalues and eigenfunctions 

In this section, we use the spectral decimation to derive finer estimates of the eigenvalues and 
the eigenfunctions, which will be used in Section]^ Of particular importance is the constant Co := 
R'{0) = the renormalization factor for the eigenvalues {Xn,k}n- Its significance derives from 
the following fact. 

Proposition 4.1. For each fc G NU {0}, the //m/t lim„^oo[-R^(0)]”^n,fc exists. 

Proof. Let p be the lowest branch of R~^, which we regard as a function on C. Via a power 
series expansion, we see that (p{z) has an attracting fixed point at z = 0 , with (y?( 0 ) = 0 and 
(/^'(O) = [i?'(0)]“^ = 2^1^ < 1- By Koenigs’ theorem (see for example |l46l §8]), the renormalized 
iterates (z i-A [(p'{0)]~'^(p'^{z)}n converge uniformly on compact subsets of a local neighborhood 
of 0 . 
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Now given a fixed k, the recurrence relation ( |3.42| ) implies that there exists no = no{k) such that 
^n+i,k = V^{^n,k) for all 71 > tiq. Combine this with the foregoing result and we conclude that the 
limit 

(4.1) lim B!{(}T\n,k = b'(0)]-”° lim [v 7 '( 0 )]-"+"V"-"°(A,„,) 

n—)-oo n—)-oo 

exists. □ 


In what follows we denote := lim„_j,oo Cq A„ fc. Our next result gives an upper and a lower 
bound on A^. 


Theorem 4.2. Fixp G (0, ^), and let k and hq be as in the proof of Proposition \4.l\ Then 

3{pq) 


(4.2) 


A 


no, A; 


1 + 


{2+pqf 


\ \ ^ ^ \ (\ + q) 

Ko,k ) < ^ A„o,fcexp - p) 


Proof As in the previous proof, let p be the lowest branch of . The lower bound on A^ will come 
from the Taylor approximation to p, while the upper bound will come from a quadratic function 
which is at least as large as p. 

Lower bound. We compute the Taylor series expansion of p about 0 to 2nd order in . 2 . It is 


(4.3) 


a(z) 


PQ KpqY 2 

2+pq^ (2 + pg)3^ ■ 


This is explained by the fact that the first derivative of the inverse function p{z) is given by p' = 
1/R\ and its second derivative is given by p" = —R"/{R!)^. Computing these derivatives at zero 
gives the quadratic function ( |4.3| ). 

Furthermore we claim that a{z) < p{z) for all .s G (0,2). It is enough to check that p'"{z) > 0. 
Here we use the identity 


(4.4) 


^ +3 (^y f ^y 

dx^ dy^ \dx) ydy"^J \dxJ 


which in our context reads 


(4.5) 


p''\z) = -R''\p{z))[p\z)r + 3[R'\p{zmp\z)f. 


Since (p'{z) > 0, we can factor out [(p\z)]^ from (4.5), and use the identity pfz) = [R'{p{z))] 
so that we reduce the original sign question to checking the sign of 


-1 


(4.6) -R'%p{z))R'{p{z)) + 3[R'\p{z))r 

(4.7) = [-6 • {3[p{z)f - Qp{z) + (2 + pq)) + 3{Qp{z) - 6 )^] 

(4.8) = ^ ■ 3 ■ {lh[p{z)f - 30p{z) + 13 - pq) 

KpqJ 

(4.9) = • 6 • [lh{p{z) - 1)^ + (1 - pq)] , 


which is always positive. This shows that p'"{z) > 0, and thus a(z) < p(z) for 2: G (0,2). 
Combined with the fact that the functions 2 i-A a{z) and 2 i-A p{z) are both monotone increasing 
on [0, 2), this implies that for each n G M, a”(z) < p"'{z) for 2 G [0, 2). 
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Fix k E NU {0} and no = no{k) as in the proof of Proposition |4.1[ Put zq = and define 
the sequence of numbers {zn}n inductively by Zn+i = (p(zn). Then, by the inequalities above, we 
have 

3(pqf 2 


(4.10) 

(4.11) 

(4.12) 


<f(Zn) > a(Zn) = 


pq 


-Zn + 




2 + pq (2 + pqY ” 

pq ^ A , 3{pq) 


1 


^ *^(^n—l) 


2+pqJ \ ' (2 + pqY 

pq \ f-, , 3{pq) 


2+pq 


1 + 


Iterating this process we arrive at the estimate 


(4.13) 


Noting that 1 


(4.14) 




no 


pq 


2+pq 


n—no n 


no 


J=no 


(2 + pg)2 

3{pq) 


*^(^n—l) 


(2 + pq) 


2 J 


3(m) 

{2^pqy 


■z > 1 , we obtain a slightly crude but still efficient estimate 


+ iZn) > Z, 


no 


pq 


2+pq 


n—no 




(2 + pq) 


2 no j 5 


which is the claimed lower bound in ( |4.2| ). 

Upper bound. To bound p{z) from above, we construct a quadratic function h{z) such that 
h{0) = 0, h{2) = p and h'{0) = (/?'(0). A simple calculation shows that 


(4.15) 


h{z) = 


pqz 
2 + pq 


4g 


and h{z) > p{z)forz E [0,2], Using this, along with the fact that = p{zn) < 

and Zn+i < \zn, we get the following estimate: 


yq 

2+pq' 


Zn \ 3 + 


p(2+g) 

4q 


Z 


(4.16) 


(4.17) 


(p{Zn) < Z, 


no 


pq 


2+pq 


n—no n 




— Xnojk 


pq 


0 =no 
n—no n 


Aq 


.‘^+pq, 

Therefore, using the inequality 1 + x <e^, 
'2 + pq 


3=no ^ 


Az, = lim 


< lim 

n^oo 


pq 

2+pq 




< 


pq 

2+pq^^° 


pq 


2+pq 


n—no n 

X, 


mo^k 




3=no 


Aq 


(4.18) 


pq 

2+pq ^ 

pq 


\ , p (2 + g) 1 ^ 

'^no,k^^P [ ^ _ p^no.k 

\ , p(2 + g) 1 , 

^no,fc6Xp[ 2 — 


This gives the claimed upper bound in ( |4.2| ). 


□ 
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For a function / : —?• E, we denote its sup norm by ||/||n,oo = sup{|/(t)| : x G Vn}. 

Likewise, the sup norm of /i : / ^ E is denoted by ||/i||oo = sup{|/i,(a:)| : x E I}. Our next result 
is an estimate on the sup norms of the eigenfunctions of A^. 

Lemma 4.3. Fix p G (0, |). Let fn^k be the eigenfunction corresponding to the {k + l)th lowest 
eigenvalue Xn^k of An - Then for every m > n > uq, 


(4.19) 


IIIItojOO a II/n,fc ||n,oo n 1 + 

j=n+l 




q — X 


n+l,fc 


In particular, if fk = limj_j,oo fno+i,k Theorem 3.3 then 

Xn,k 


(4.20) ||//c||cx) — II ./n,/c ||n,c» 6Np 
Proof Let us introduce the function 

(4.21) F{A,B,z) = 


P 


Q Xfi-\-i^k 2 p 

q{l — z)A + pqB 
{q-z){l+p- zf 


which is derived from the eigenfunction extension algorithm ( |3.1| ) and (3.2). Note that .2 G [0,p) 
because the extension uses the lowest branch of R~^ starting from level hq. First we would like to 
control the linear growth of z ^ F{A, B,z) — F{A, B,0): 

\F{A,B,z)-F{A,B,0)\ 

<\A\ 


1 

1 

+1^1 

pq 

pq 

1 

+ 

1 

1 

to 

{q - z){l + p - z) 

1 _ p 2 


<max(|A|, \B\)q 
(4.22) 

= max(|^|, \B\)q\z 


(1 — z)(l — p^) — 1 + 2z + 

+ p 

2z — z^ 

{q - z)(l +p- z)(l -p2) 

(q- 

+ 

1 

1 

to 


1 + p"^ — z 


{q — z){l +p — z){l - p^) 


p 


2-z 


{q - z){l + p — z){l - p^) 


Since 0 < z < p < |, the absolute value terms in the RHS of (4.22) are positive, so we can drop 
the absolute value signs and add the two terms in the bracket to get 


\F{A,B,z)-F{A,B,D)\ 

^ n/II I DU I 1/ l + 2p + p2-z-pz 

< max(|A|, |i5|)g|z| 


(4.23) 


= max(|A|, \B\)q\z\{l F p) 
= max(|^|, \B\) 


(g-z)(l+p-z)(l-p2) 

1 + p — z 


{q — z){l + p — z){l - p^) 


{q-zf 


which implies that 

(4.24) |F(A5,z)| <max(|Al|,|S|)|z|^^ + |F(A5,0)|. 

[q-z) 
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We use ( |4.24| ) to estimate the sup norms of the eigenfunctions: for all n > no, 

ll/n+l ,A:||n+l,oo ^ ||/n,fc||n,oo 

^ ||/ri,fc||n,c 


l^n+1I 


(4.25) 


^ \\fn,k\ 


Q ^n+l 
^n+1 


+ |F(1.1.0)1 

+ 1 


^n+1 

Iterating the inequality (4.25) and using the fact that Zn+m < (f )"* (3.14) gives 


II/tojA: ||m,oo ^ ||/n,fc|U,oo f 1 

j=n+l 


Zn 


,,k ||n,oo 


(4.26) 


< ll/n,, 

||/n,/c ||n,cx) 


n 1 


j=n+l 

m 

n 

j=n+l 


<1 - 

p\j-^ Z„ 


q - Zn+l 


^ I PW ” ^n,k 

22 Q /j 


for all m > n > no. This shows ( |4.19| ). 

Next, using the triangle inequality and taking the supremum, we have 


(4.27) 


sup|/fc(a:)| - sup |/^,fc(a:)| 


X^I 


x£l 


< sup|/fc(x) - fm,ki^)\- 


X^I 


Recall from Theorem |3. 3 [that the limit limj^oo fno+i,k is uniform on /. This along with the bound 
( |4.27[ ) implies that 

(4.28) ll/fclloo lim II/rn,fc ||m,oo 

m^oo 

So by taking the limit m ^ oo on both sides of ( |4.3| ), we arrive at the estimate 

' P\^ ^ ^n,k 


Wfi 


IIoo 


< U 


n,k ||n,oo 


lim IT ( 1 + 

m. — VcxT) \ 


q — X 


(4.29) 


< ||/n,fc||n,ooexp j {^- 

\j=n+l 
^n.k 


j=n-\-l 

oo . A 

P\^ ^ Xn,k 
22 q /j 


n+l,/c 


||/n,fc ||n,c» *3Xp 


P 


□ 


A„+i,fc2-py ’ 

where in the second line we used the inequality 1 + a: < e^. This proves ( |4.20[ ). 

The following result provides a quantitative estimate of the convergence of fn,k to fk in sup norm. 
As in the previous section, we harmonically extend fn,k from 14 to I, and abuse notation by calling 
the extension fn,k still. Then A^fn^ki^) = 0 for all a: G / \ 14- 

Theorem 4.4. Fix p G (0, |). Then 


Wfk /njfcllcc ^ Cq A^;II14^00 ||g||ooexp 


A 


n,k 


p 


q - An+l,fc 2-pJ ' 


(4.30) 
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where g \ I x I ^ is the Green's function associated with and 

(4.31) Ikiloo := sup g{x,y). 

{x,y)GlxI 

Proof. Since fk and fn,k agree on 14, it is enough to estimate their difference on / \ 14- Based 
on the construction described in Section]^ / \ 14 is the disjoint union of {F^„{(0,1)} : \w\ = n}, 
where w = wiW 2 • • • is a word of length n with w* e {1,2,3} for 1 < i < n, and Fy^ := 
Fwi ° o • • • o Fyj^. So our task is to show that for every word w of length n, the function 
\{fk — fn,k) o on (0,1) has a uniform upper bound. 

Our strategy is to exploit the self-similarity of the fractal Laplacian (Proposition |2 .4| ) , as well 
as properties of the corresponding Green’s function. We remind the reader that is the 

Green’s operator associated to A^. It admits an integral kernel g{-, •) called the Green’s function, 
defined by 

(4.32) {Gt,u)iy) = g{y, y')u{y') dy{y'). 

The existence of and g follows from the theory of Kigami llTTl §3.5-§3.6]. In particular, G^ : 
C{I, /i) C{I), and p is a nonnegative continuous function on / x /. 

To begin the proof, we start with the self-similarity of A^ (Proposition |2.4[ ): 

(4.33) A^((/fc - U,k) o Fy,) = CoA^(/fc - /,,,) o Fy, on /. 

Since fn,k is harmonic with respect to A^ on F^„{(0,1)} for every word w of length n, it follows 
that 


(4.34) ^u,ifk-fn,k)°Fy, = Af,fkoFy, = XkfkoFyj on (0,1). 
Combine ( |4.33| ) and ( |4.34| ) and we get 

(4.35) A,{{fk - fn,k) o Fy,) = Co^Xkfk o Fy, on (0,1). 

Now apply the Green’s operator G^ on both sides of ( |4.35[ ) to get 

(4.36) ifk - fn,k) o Fy, = C^^XkG^ifk o Fy,) on (0,1). 
Using the representation ( |4.32| ) we obtain the estimate 

\ifk - fn,k) o Fy,\ < Co”Afc||G'^(/fc O F^)||oo 


(4.37) 

This proves that 

(4.38) 


<Co^h\\g\U\ifkoFy,)\\^ 

< Co”Afc||5f||oo||/fc||oo on (0,1). 


Wfk - /n,fc||oo < Co”Afc||5f||oo||/fc| 


In order to turn the RHS of ( |4.38| ) into a mo re use ful estimate, we use (4.20) to replace ||/; 
by ||/n,fc||n,ooexp which yields (4.30). 


II oo 
□ 
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5. Estimates on the solution oe the wave equation 

We now apply the results of Sections|^and|^to estimate the solution of the wave equation on the 
interval / endowed with the fractal measure fi. 

Numerically, we can only compute the eigensolutions of the fractal Laplacian up to a finite level, 
so in practice we solve the “approximate” wave equation 

= onK, 

dtUn{-,0) = 0 onVn, 

3^0 

where 0'kfn,k is the approximate (5-function built up from the first \Vng \ = (3**° -|- 1) 

k=0 

eigenfunctions of A„ (with Anfn,k = ^n,kfn,k)^ and ak := Ono.fc A 0 are the coefiicients found in 
Section |33| Throughout the section no will be fixed, and we will not mention no explicitly unless 
the context demands it. 

Following the exact same argument, the solution to ( |5.1| ) has the series representation 

3'*o 

(5.2) Un{x, t) = oikfn,k{^) COS for all x G 14 and f G [0, T]. 

/c=0 

For each t, we harmonically extend the function x i-A Un{x, t) from 14 to /. This procedure allows 
us to compare Un{x, t) with 

3"o 

(5.3) u{x, t) = Oikfkix) cos (^t'i/Xk^ for all a: G / and t G [0, T], 

k=0 

the solution of the wave equation on (/, n) whose initial condition is the truncated series repre¬ 
sentation of the (5-impulse. We note that u is differentiable in t and continuous in x because the 
eigenfunctions functions fk are continuous. However it is highly localized function at t = 0, and 
therefore it mimics wave propagation from a delta function initial values. 


5.1. Convergence of approximate solutions of wave equation. In this subsection we establish an 
upper bound on 


(5.4) 



— u 




for all n and t, uniform in x. This would then give us the convergence of Un to u at fixed t and 
uniformly in x. Note that we are normalizing t in such a way that the orthogonal projection of the 
wave onto the lowest eigenfunction (corresponding to eigenvalue o) propagates at speed 1. 


Theorem 5.1. FixpE (0,|). Let Un and u be respectively defined as in (5.2) and (5.3). Thenthere 
exists a positive constant C = C{no,p) such that for each t G [0, T] and n > uq, 


Un 



— U 



til 


\/Ai 


(5.5) 


sup 

xel 


< C (f V 1) Co 
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Proof. Using the series representations ( |5.2| ) and ( |5.3| ) and the triangle inequality, we find that ( |5.4| ) 
is bounded from above by 


(5.6) 


f \A 

fn,kix)cOS tlT^ 


k=0 




- fk{x) cos ( tir 


7^ 


which, by a simple manipulation using the sum-to-product trigonometric rules, is equal to 


S’^o 


(5.7) 


I - fk{x)] COS cos 

+ fk{x)] sin {tTTA+^ sin (tTrA” J | , 


k=0 


where A^^ := ^ i \/^) ‘ again the triangle inequality, we can estimate (5.7) from 

above by Ii + I 2 , where 


3'*o 


(5.8) 


(5.9) 


k ■= \fn,k{^) - fk{x) \ • |cos (fTrAjfc) I |cos (f7rA„ fc) | , 

fc =0 

3'*o 

h ■- \fn,k{^) + fk{x)\ • I sin (fvrAjfc) | | sin (fvrAj^ I . 


k=0 


The key term to control in Ii is | fn^k (^) ~ /fc (^) I > while in I 2 it is | sin {tirA^ j.) |. For the former we 
invoke Theorem |4.4[ while for the latter we apply the Taylor expansion 


(5.10) 


sm 


{tT^K,k) I < |^n,fc| + O {{tAn,k?) . 


For terms other than these two, we apply the simple minded estimates | cos I < 1 , |sin (tTTAjfc) I < 

1 , and 

\fn,k{^) + fk{x)\ < 2|/„,fc(a:)| + \\fn,k - fkWoo 

(5.11) <2||/„,fc||oo + 0(Co”). 

First let us estimate A”^., which amounts to controlling the ratio ^ 7 . By Theorem 

1 + ^ ^ CXp (Ah,!^^) 


4.2 


(5.12) 


(2+pq)^ ' 


exnlA 

exp \An,k 2q{2-p) 


A, 


n,l 


. — < 
Afc 




(2+pq) 


2 '^n,k 


We note that in this case we have 0 ^ k ^3'^° < n and so Theorem [43| is applicable. 

In what follows we denote Do := 2 q( 2 -p) Observe from the discussion in the 


(2+pg)2 ' 


proof of Theorem 4.2 that Do > Di whenever p G (0,1]. 

From Proposition 4.1, we know t hat fo r each fixed k, A^ ^ = ( 9 ( 00 **) as n —?• cx). Thus upon 
expanding the LHS and the RHS of (5.12) up to the 0{Cq^) terms, we get 


1 + 2 ~ DoA„^fc) +o(Co”) < 


^n,k Ai 
An,l Afc 


— 1+2 (^oAn.i “ DiA,i^fc) + o(Cq**). 


( 5 . 13 ) 
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It then follows that 


(5.14) - (DiA„,i - DoA„,,) + o(Co") < K,k 

Plugging this into ( |5.10| ) yields 
I sin (tTrA-^) I 

- ^t/^niax(|DoA„,i - DiA„,fc| 
4 V 


^ < 7 (DoA„,i - DiA^,,) + o(Co”). 


|DoA„^/c — DiA„^i|) + ©(tCo"") 


(5.15) 


= ^./5|DoA„,fc-DiA„,i|+o(tC 

4 V Ax 


-n\ 
0 ) 


Putting all the estimates into (5.8) and (5.9) gives 


3"o / A 2 A 

h < Co”|| 5 f||ooy^«fcAfc||/„,fc||„,ooexp (-—--) + o (Cq”) . 

^ \q-K+i,k^-pJ 


(5.16) /2<y5]«fc||/, 




k=0 


^|DoA„,fc-DiA„,i| + o(tCo”). 

Al 


This means that (5.4) is bounded above by 


Sno 




n.K oo 


k=0 


AfcCo^llfi'llooexp 


A 


n,k 


(5.17) 


tlT / I .. .. I 

~~2\ Y “ DiA,i^i| 


q - \n+i,k 2 - p 

+ o ((t V 1) Co”) . 


It remains to explain how ( |5.5| ) follows from ( |5.17| ). For 0 < A: < 3**°, A^ < A 3110 . Also by 
Proposition |4.1[ there exists a positive constant C independent of n such that 

(5.18) |DoA„,fc - DiA„,i| < CCriDoA, - DiAi|. 

By an argument of Kigami OtII . the Green’s function corresponding to on I can be con¬ 
structed independently of the measure p, whence independently of the value p. In particular when 
p — we recover the classical Green’s function on I with Lebesgue measure, g{x,y) = (a: A 
y){{l - x) A (1 - y)). Thus ||p||oo = 

For the exponential in the first term, note that since p G (0, |), 


(5.19) 


A 


n.k 


< 


A, 


n.k 


4 A^ 


n,k 


Q 2 p 2 ^n-\-l^k 


2—2 3 2 — An+l,fc 


By Proposition |4.l[ the RHS of(5.19)is(9(CQ”) asn —?• 00 . Thus the exponential is exp(0(Co ”)) = 
1 + 0(C„-"). 

Finally, we claim that there exists a constant C such that 


S^o 

^ ) ^k II ./n,fc II00 A C 
k=0 

for all n > hq. To explain this, note that = T by virtue of our choice of the initial 

eigenfunctions (/o,i and /o, 2 ) and initial weights aop = 00,2 = and the matching condition 
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Figure 5 . 1 . A graph of y{x) = otk \ fn , k { x )\ and in order: uniform spacing of 
points, uniform resistance between points, and uniform measure. 


( |3.33[ ). Then the rest of the proof follows from Section 8 of L. Rogers’ paper IBTII . Numerically, C 
is slightly above 1 , see Figure [5?Tj □ 


6. Numerical computation oe eigenfunctions and solutions of the wave equation 

We present some of the numerical results obtained by our spectral decimation method. The 
spectral decimation is an iterative method, and the code repeats the calculation done in section 
This code, which is used to produce pictures and to perform the experiments, and a graphical 
user interface to recreate the results can be found at http: //homepages. uconn. edu/f ractals/ 
fractalwave/, Here we give a representative variety of figures detailing some of the numerical 
calculations that have been performed. Figures |TTj and |3.2| show the first 25 eigenfunctions, and, in 
particular, the ways in which the symmetry is broken for values of p and 1 — p. Figures and |3.4| 
show the quality of the approximation for the delta function for various values of rig. In particular 
one can see that for small values of p (in our case p = .2) the approximation is significantly better 
than for the corresponding values of 1 — p (in our case q = . 8 ). This shows why our efforts focused 
on the cases where p < .5. The next set of figures ( |6.1[ |6.2[ |6.3[ |6.4| ) heuristically suggest that 
the visible portion of the wave propagates at a speed proportional to . But, further investigation 
will be needed to show this more precisely. Figures [63] shows three different parametrizations of a 
representative eigenfunction. 
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Figure 6 . 1 . Uniform spacing of points at time t 
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Figure 6 . 2 . Uniform spacing of points at time t 
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Figure 6 . 3 . Uniform spacing of points and, from left to right: time t = 0.1, 0.2, 0.3, 
0.4. with p = 0.1 
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Figure 6 . 4 . Time t = 0.4 and p — 0.1 and from left to right: uniform spacing of 
points, uniform resistance between points, and uniform measure. 
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Figure 6 . 5 . Eigenfunction 53, withp = 0.1 and in order: uniform spacing of points, 
uniform resistance between points, and uniform measure. 
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